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ABSTRACT 


-Ibis report describes findings of research conducted at the University of 
Nebraska Agricultural Research and Development Center in 1985 and during a 

pilot study at the FIFE site in 1986. 

The primary objectives of the 1985 study were (1) to test the feasibility 
of using radio frequency receivers to collect data from automated weather sta- 
tions and (2) to evaluate the use of the data collected by the automated 
weather stations for node ling the fluxes of latent heat, sensible heat and 
radiation over wheat. The model Cupid was used to calculate these fluxes 
which were compared with fluxes of these entities measured using raicrometeoro- 

logical techniques. 

The objectives of the 1986 study were (1) to measure and model reflected 
and onitted radiation streams at a few locations within the FIFE site and (2) 
to compare nodeled and Treasured latent heat and sensible heat fluxes frcm the 

prairie vegetation. 

During the 1985 study we found that it was feasible to use radio fre- 
quency transmitters and receivers to collect data frcm automated weather sta- 
tions. We found that nodeled estimates of latent heat (LE), sensible heat (H) 
and net radiation (Bn) agreed well with the fluxes measured over wheat, espe- 
cially during the time when the wheat crop was growing vigorously. The 
agreement between treasured fluxes and modeled estimates diminished later in 
the season as wheat an thesis ended and the wheat began to senesce. 

Reflectances and transmittances from leaves of some of the more abundant 
grass and forb species were made using a specially designed radicmeter. This 
instrument was equipped with filters to measure values in the same waveband 
intervals as those found on the Barnes Modular Multiband Radiometer (MMR) . 
All grass and forb leaves exhibited peak reflectance and transmittance values 
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in the visible region in SWR band 2. Peak reflectance in the near infrared 
region occurred in MIR band 4 for grasses, but for fort) leaves peak reflec- 
tance occurred in M4R bands 4 or 5. The magnitudes of reflectances and 
transmittanoes were similar for measurements made on either the adaxial (top) 
or the abaxial (bottom) surfaces of grass leaves, but such was not the cas e 
for forb leaves. These results suggest that it should be possible to charac- 
terize the leaf reflectance and transmittance by making measurements on either 
the adaxial or abaxial surface of grass leaves, but transmittance and reflec- 
tance measurements must be made on both sides of forb leaves. 

Reflectances from the prairie vegetation canopies were measured in the 
fWR wavebands at two sites, one on the Konza (a sloping site) and the second 
off the Konza (a level site). Reflectances were made at seven angles in the 
principal plane of the sun. The angles were -50, -30, -20, nadir, 20 , 30 and 
50. Negative values signify backscattered radiation and positive values for- 
ward scattering. There were several similarities in the reflectance patterns 
made at different locations, but also some differences. At both sites, the 
minimum reflectances tended to occur at 20 or 30°. On the Konza site the 
maximum reflectance occurred consistently for MMR band 4, but at the off -Konza 
site, it occurred for band 4 or 5, depending on the specific sampling loca- 
tion. Forward scattered reflectances were generally 2-5 percent lower in 
absolute terms than backseat ter reflectances. The effect of changing zenith 
angle on reflectances was inconclusive over the range of solar zeniths we 
experienced. Changes in canopy or leaf characteristics during the measurement 
periods may have masked any solar zenith effects. 

Global albedo values were calculated from bidirectional reflectances in 
the MIR wavebands using weighting coefficients determined from solar radiation 
data by Moon (1940) and hemispherical reflectances calculated by the model of 
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Walthall (1985). Calculated values ranged from about 20-25%. These values 
are reasonable, but albedo measurements using total hemispherical albedcroeters 
should be made concurrently with the bidirectional reflectance measurements 
for comparison purposes. 

Measurements of fluxes of IE, H and soil heat (S) were made on eleven 
days over wheat in 1985 and two days over the grasslands in 1986 using Bowen 
ratio energy balance and eddy correlation techniques. The partitioning of net 
radiation into sensible heat was dependent on soil moisture, atmospheric con- 
ditions and the stage of vegetative development. The measured values of these 
fluxes were used to compare with fluxes estimated with the model Cupid. 

The performance of the Cupid model for estimating l£, Rn and H fluxes in 
1986 from prairie vegetation was poorer than that for the more homogeneous 
wheat canopy in 1985 even at comparable LAI values. Modifications to Cupid 
are needed to account for the nonhcmogeneous, randomly-spaced clumps of 
prairie vegetation as well as a need to incorporate the influence of leaf 
litter on the fluxes of IE, H and S. 
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introduct io n 


The main objective of the International land Surface Climatology Project 
( ISLSCP) has been stated as "the development of techniques that may be applied 
to satellite observations of the radiation reflected and emitted from the 
Earth to yield quantitative information concerning land surface climatological 
conditions." To accomplish this objective a major field study called Fire - 
the First ISLSCP Field Experiment - was conducted in 1987 on and near the 
Konza Prairie in Kansas. Daring that experiment two types of measurements 
were made, (1) long-term monitoring of basic meteorological parameters and (2) 
extensive surface and aircraft measurements made during four intensive field 
campaigns (IFC's). Among activities taking place during these IFC's is the 
determination of radiation fluxes, sensible and latent heat fluxes and the 
measurement of selected biophysical properties. This document describes 
results from pilot studies conducted in 1985 at the University of Nebraska 
Agriculture Research and Development Center and in 1986 at the FIFE site in 
Kansas to collect information and develop techniques contributing to success- 
ful measurement of selected parameters from these two measurement types. 

Oar basic overall objective in conducting these pilot studies was to 
document and understand certain important physical and biological conditions 
and processes which occur on the prairie site. Tb do this we made selected 
meteorological, plant and soil measurements to permit us to compare flux esti- 
mates with a plant simulation model, Cupid, with fluxes of radiation, sensible 
heat and latent heat calculated using micrometeorological techniques. Speci- 
fically, our objectives were (1) to determine the feasibility of using radio 
frequency transmitters and receivers to collect data frcm automated weather 
stations that will aid in determination of the me so- and micro-climatic 
variability at the FIFE site; (2) to measure and model reflected and emitted 
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radiation streams from various areas within the prairie site; (3) to make 
reliable estimates of sensible and latent heat fluxes over the prairie 
employing the eddy correlation technique and other micrcmeteorological methods 
and then to (4) evaluate the model Cupid for estimating the radiation, sen- 
sible and latent heat fluxes from routine atmospheric, plant and soil measure- 
ments* Such models offer a means to provide areal estimates of these fluxes 
from a limited number of point measurements. 

Procedures and results of our research efforts on each of the four objec- 
tives specified above are presented in this report. 

Objective (1) - Automated Data Collection Through Radio Frequency 
Cownunication “ 

In 1985 equipment was procured to perform automated data collection via 
Radio Frequency (RF) communication. Daring the summer of 1985 the equipment 
was installed at the Agricultural Meteorological Laboratory at Mead, Nebraska 
for testing. Testing was carried out with an on-site IBM computer and an 
existing automated weather station at a nearby, permanent location. Linkage 
between the computer-modem-RF radio and the data logge r-modem-RF radio was 
established and tests indicated that communication was satisfactory. 

In 1986 additional equipment was purchased so that the field sites could 
be located on the Konza prairie. Three stations were installed to collect 
data in late July and early August 1986. A computer was installed at the 
Konza Prairie Headquarters building and an RF antenna was deployed on the 
roof. Descriptions for the three weather stations are given in Table 1. 
Stations #1 and #2 were located on the Konza Prairie while #3 was located just 
southeast of the Konza. Although the computer was only checked periodically 
(one to two times per week), the data collection process was reasonably suc- 
cessful as indicated in Table 2. The Konza #1 station utilized a CR21X data 
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Table 1. Description of the 3 surface weather stations deployed in 1986 


Instrument 

Konza #1 
426 

Konza #2 
400 

Konza #3 
400 

Data logger 

CR21X 

CR21 

CR21 

RF-modem 

yes 

yes 

no 

RF-rad io/an tenna 

yes 

yes 

no 

Tape recorder 

no 

no 

yes 

Anemometer 

yes 

yes 

yes 

Wind vane 

yes 

yes 

yes 

Pyrancmeter 

yes 

yes 

yes 

Air temperature sensor 

yes 

yes 

yes 

Air relative humidity sensor 

yes 

yes 

yes 

Soil temperature probe 

yes 

yes 

yes 

Precipitation gauge 

yes 

yes 

yes 

Net radiometer 

yes 

no 

no 

PAR sensor 

yes 

no 

no 

Canopy temperature sensor 
- Nadir view 

yes 

no 

no 

Canopy temperature sensor 
- 30° view angle 

yes 

no 

no 


- 6 - 






Table 2. Time period of data collection. 


7/21 


7/28 8/4 8 / 11 


Konza #1 


<- 


->l 


Konza #2 


|< >1 1 < >1 >! 


Konza #3 


<- 


•>i 


Experience with the CR21 in unattended rode would indicate that it worked 
well ^ recorder (.3) but not with the EF-linkage <»2>. the CR21X 
(#1 above) and RF-linkage gave an excellent performance. 
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logger which performed extremely veil in the RF-linkage. Konza #2 utilized a 
CR21 data logger and would not operate continuously in RF-linkage. Konza #3 
also had a CR21 data logger but was not part of the radio ocranunication test 
(data were recorded on a cassette tape). The CR21X data logger manufactured 
by Campbell Scientific of Logan, Utah was conveniently employed, with no modi- 
fications, to collect sophisticated surface weather measurements on the Konza 
Prairie. We concluded that the CR21X data logger could be easily interfaced 
to additional sensors and used with RF-ccmmuni cat ions to establish a reliable 
weather network for FIFE. A CR10 data logger with RF linkage is now avail- 
able. It is less expensive than the CR21X and could be used to obtain the 
parameters measured in this study. 

Objective (2) - Measuring and Modeling Reflected and Bnltted Radiation 

Data collected to help meet this objective were obtained during late July 
and early August of 1986 at the FIFE site. We designed the study, in part, to 
help familiarize us with the topography and vegetative composition of the site 
and to test and develop equipment and procedures to be employed during the 
1987 FIFE study. The major focus of our effort was to obtain reflectance and 
transmittance data fron sane of the dominant vegetative species and ooimuni- 
ties located within the FIFE site. 

Instrumentation used to gather data on reflectance and transmittance 
included the Barnes model 12-1000 modular nultiband radiometer (MMR) , and a 
LI-OOR model 1800-12 integrating sphere used with a specially designed radio- 
meter. The Barnes mR was used to collect reflectance data on vegetative com- 
munities along transects that may have differences due to topography- 
chief ly, slope and aspect. The integrating sphere was employed to collect 
data on reflectance and transmittance for individual blades of grass or leaves 
of forbs. Measurements of leaf area index (LAI) by destructive and nondes- 
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tructive methods along with sampling for vegetative dry matter were done for 
some of the locations where radiation measurements were made. 


A. Leaf Transmittance and Reflectance : 

Modeling of radiation fluxes requires knowledge of the transmissivities 
and reflectivities of leaves of the dominant plant species growing on the FIFE 
site. We measured the reflectivity and transmissivity of the leaves of 
several grass and forb species using a LI-COR integrating sphere fitted with a 
specially designed radiometer to collect data in the wavebands of the Barnes 
MMR. These wavebands are: (1) 0.45-0.52 pm, (2) 0.52-0.60 pm, (3) 0.63 0.69 
pm, (4) 0.760.90 i*n, (5) 1.15-1.30 pm, (6) 1.55-1.75 pm, and (7) 2.08-2.35 pm. 
Leaves of forbs or blades of grass were placed into the integrating sphere 
without detaching them from the plant; that is, the plants were studied 
nondestructively . 

Of the hundreds of species of grasses and forbs found at the FIFE site, 
eleven of the rrost oommon ones were studied (six grass species and five forb 
species). The species name, common name of that species, whether it is a 
grass or forb, and the number of samples of that species for which data were 
collected are listed in Table 3. 

Reflectances from the adaxial and abaxial (rft and rfb, respectively) of 
a blade of Andropogon qerardii as well as the transmittances through the top 
and bottom of the blade (tft and tfb, respectively, are depicted in Figure 1). 
The transmittances through both the top and bottom of the blade were essen- 
tially the same and were lower than the reflectances for the first three wave- 
bands. In band 4 reflectances and transmittances were very similar for either 
side of a leaf. In bands 5 and 7 the reflectances from adaxial or abaxial 
surfaces tended to be within one or two percent (absolute values) of each 
other and were generally lower than the transmittances. Transmittance through 
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Table 3. Vegetation studied at the FIFE site 


Species 

Ccranon name 

Species 

Type 

Number 

or 

Samples 


Andropoqon gerardii Vitman 

Big blues tem 

Grass 

3 


Panicum virgatum L. 

Switchgrass 

Grass 

2 


Agropyron smithii Rydt. 

Western wheatgrass 

Grass 

2 


Eragrostis spectabilis (Pursh) Gread. 

Purple lovegrass 

Grass 

1 


Bouteloua curtipendula (Michx.) Ibrr. 

Sideoats grama 

Grass 

1 


Sorghastrum nutans (L. ) Nash 

Indiangrass 

Grass 

1 


Salvia pitcheri Tbrr. 

Pitchers sage 

Forb 

1 


Afemonia baldwinii Ibrr. 

Iconweed 

Forb 

1 


Ambrosia psilostachya DC. 

Western ragweed 

Forb 

1 


Ceanothus americanus L. 

New Jersey tea 

Forb 

1 

- 

Rhus glabra L. 

Sumac 

Shrub 

1 
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Figure 1. 


Waveband Number 

Leaf reflectance from the adaxial (rft) and abaxial (rfb) and tranamittance through the adaxlal (tft) and 
abaxial (tfb) in the MMR wavebands for Andropogon gerardii. 





the top of the blade was essentially the same as reflectance from either the 
top or bottom of the leaf for band 6* Overall 9 the reflectances and transmit** 
tances in bands 1 through 3^ which fall within that region of the electromagne** 
tic spectrum known as the visible or photosynthetically active radiation 

A 

(PAR) are nuch lower than those found in the remaining wavebands. In the PAR 
region actively growing plants are highly absorptive. Fteflectances and 
transmittances in bands 4 through 7 fall within the infrared region and are 
much higher than those in the PAR region. 

Three samples of A. gerardii exhibited similar reflectance character is- 

l 

tics in all bands (Fig. 2). In the PAR region reflectance and transmittance 
values peaked in band 2 (the green region), while in the infrared region 
reflectances peaked in band 4 (0.76-0.90 pm) and transmittances peaked in band 
^ (1.15—1.30 iin) . This pattern was true for all grasses. 

Panicum virgatum (Fig. 3) exhibited trends very similar to A. gerardii . 
However, P. virgatum transmittances were slightly higher in all bands, except 
band 3, than A. gerardii . Reflectances for the two samples of P. virgatum 
behaved similarly, with one leaf showing consistently higher reflectance in 
all wavebands than the other (Fig. 4). 

Data for Agropyron smithii (Fig. 5) depict values for transmittance and 
reflectance that were higher than those for P. virgatum but lower than A. 
gerardii in wavebands 1 through 3. In band 4, A. smithii peaked at higher 
reflectance than any other grass evaluated during the pilot study. Reflec- 
tances in bands 5 through 7 for A. smithii were greater than those observed 
for the previous two species. Reflectances for the two samples of A. smithii 
were very similar (Fig. 6). 

Values of reflectance and transmittance for Bouteloua curtipendula are 
plotted in Fig. 7. Transmittances were much higher than reflectances in the 
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REFLECTANCES FOR TWO SAMPLES OF 
PANICUM VIRGATUM 



FTggre 4. Comparison of leaf reflectances for leaves from 2 Panicum virgatum plants. 
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Comparison of leaf reflectances for leaves from 2 Agropyron smithii plants. 
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figure 7. As in Fig. 1 for Bouteloua curtipendula. 


longer wavelengths, a pattern also observed for Eragrostis spectabilis (Fig. 
8). B. curtipendula and E. spectabilis had the highest transmittances in the 
longer wavelengths of all the grasses studied. Also E. spectabili s had the 
most irregular data (Fig. 8). It is not known if this is characteristic of 
this species or if the irregularities were due to measurement problems. 

Figure 9 is a plot of the data for Sorghastrum nutans. The reflectance 
in band 4 was the second highest while transmittances in the visible wavebands 

were the lowest of the grasses evaluated. 

Reflectances from the adaxial surfaces of all grass species followed the 
same general pattern; that is, peak reflectances in the visible portion of the 
spectrum in band 2 and peak reflectance in the NIR spectrum in band 4 (Fig. 
10). Seme species differ from each other in some wavebands, yet are similar 
in others. Because of this, it is important to identify the dominant grass 
species and characterize their reflectance and transmittance properties in 
order to better model the surface radiation fluxes. Our data suggest that 
measurements from both the top and bottom of grass leaves may not be required 
since transmittances and reflectances on either side are within three to five 
percent in absolute terms. 

As was the case for grasses, forbs displayed a peak reflectance and 
transmittance value in band 2 of the visible region (Figs. 11-15). However, 
in the infrared region transmittance peaked in either band 5 or band 4, while 
the reflectance usually peaked in band 4 but tended to be "flat" between bands 

4 and 5. 

In contrast to the grasses, reflectances and transmittances of seme forbs 
varied according to the leaf surface. For example, transmittances through the 
top and bottom surfaces of Vemonia baldwinii (Fig. 12) and the shrub Rhus 
qlabra (Fig. 14) were very different in the visible region but similar in the 
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Figure 8* As in Fig. 1 for Eragrostis soectabilis. 
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Figure 9. As in Fig. 1 for Sorghastrum nutans . 
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Rgure 10. Comparison of reflectances in the MMR wavebands for all grasses measured in the study. 




Reflectance ar Transmittance (%) 


REFLECTANCE AND TRANSMITTANCE FOR 
SALVIA PITCH ERI 



Figure 11 . As in Fig. 1 for Salvia pitcheri. 
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Figure 12 As in Fig. 1 for Vernonia baldwinii. 
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Figure 13. As in Fig. 1 for Ambrosia psilostachya , 
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Figure 18. As in Fig. 1 for Ceanothus americanus . 


NIR region. Reflectance from the adaxial surface of the leaf was much lower 
than from the abaxial in wavebands 4 through 6. However, seme species of 
forbs did not display such large differences in reflectances and transmittan- 
ces for the different leaf surfaces as shown, for example, by Ceanothus ameri- 
canus (Fig. 15). Reflectance and transmittance properties of the forbs, 
except for Salvia pitcher i , are very similar in the infrared wavebands, but 
are not similar to each other in the visible wavebands (Fig. 16). Forbs may 
not make up a very significant portion of a satellite pixel, but the contribu- 
tion of forb reflectance to the overall surface reflectance at various loca- 
tions within the FIFE site needs evaluation. 

B. Bidirectional Reflectance and Quittance : 

Reflected spectral radiation from prairie surfaces was measured with a 
Barnes fWR model 12-1000. It was necessary to take M® measurements using a 
portable mast because of restrictions placed on driving vehicles on the Konza 
prairie. The mast allowed us to take measurements along predetermined tran- 
sects, and enabled us to record changes in spectral properties of vegetative 
communities that may have resulted from changes in slope and aspect, or 
changes in species composition. 

The portable mast (Fig. 17) was designed to hold the PWR approximately 
3.1 m above the soil surface. At this height the radiometer viewed a spot 
size of about 0.75 m in diameter at nadir. Measurements were trade from seven 
different view angles. These angles were nadir and 20°, 30°, and 50° to 
either side of nadir. The mast was oriented so that readings were made in or 
near the principal plane of the sun. 

Data were collected on a total of twenty-three plots on transects at two 
different sites. The first site was located on the Konza Prairie southeast of 
the headquarters building in the NE 1/4 of section 24 and SE 1/4 of section 23 
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Figure 16. 


As in Fig. 10 for forbs. 





in the Swede Creek quadrangle. It is referred to as Konza 1 on the figures. 
The second transect was located on a gently sloping area near the flux site 
referred to as site 11, 16, 18 in the 1987 FIFE study. It is referred to as 
Konza 2 on the figures. These transects ware chosen to exhibit differences in 
topography and species composition. Data were taken from Konza 1 and Konza 2 
on July 29 and July 30, 1986, respectively. 

MMR voltages were recorded on an Qnnidata polycorder and converted to 
reflectance factors (RF's) using the procedure described by Gardner (1983). 
These RF values were plotted and examined as a function of view angle and also 
used in the calculation of albedoes and vegetation indices. 

1. Reflectance as Function of View Angle. 

a. Konza 1 Site : The influence of slope, time of day and species com- 
position on bidirectional reflectance factor distributions for grassland sur- 
faces in the principal plane of the sun are shown in Figures 18-23. The dips 
in the curves in the backscatter direction are probably caused by instrument 
shadowing of the surface viewed by the MMR. Without exception, backscatter 
reflectance was greater than forward scattered reflectance. The minimum 
reflectance occurred not at nadir, but rather at view angles of about 20° or 
30° for all measurements made on the north-facing slope (Figs. 18, 19, 20 and 
22) and at nadir for the south facing slope (Fig. 23). 

Reflectance decreased as the solar zenith angle decreased, especially as 
the view angles approached nadir. The effects were most pronounced for NIR 
reflectances (bands 4, 5) although visible reflectances (bands 1, 2, 3) also 
decreased with decreasing solar zenith angle and view angles near nadir. 

Some general conclusions for reflectance curves from the Konza 1 transect 
are summarized as follows: 

1) The general order of reflectance was, in order of increasing values, 
bands 1, 3, 7, 2, 6, 5 and 4. 
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figure lo. Percent reflectance as a function of view angle for plots 1 and 2 at the Konza 1 site 
on day 210 (July 29, 1986). The solar zenith angle was 41°, leaf area index was abou- 
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Figure 23. Percent reflectance as a function of view angle for plot 13 at the Konza site 1 on day 210. 

The solar zenith angle was 28°, leaf area index was about 1.5 and a south-facing slope of 15°. 



2) Reflectance magnitudes and patterns were very similar for the visible 
bands of 1 and 3 which are located in the chlorophyll absorption 
bands and for band 2 (visible) and band 7 (mid— IR). 

3) Minimum near-IR reflectance occurred at 20-30° in the forward scatter 
direction and maximum at 50° in the backscatter direction. Forward 
scatter reflectance at 50° was generally 10-20% (15-30% relative dif- 
ference) lower than backscatter at 50°. Minimum values were 15-30% 
(30-50% relative) lower than the maximum values. 

4) Visible reflectance generally decreased from the backscatter direc- 
tion at -50 to the forward scatter at 50°. Forward scatter was 
generally 2-5% lower (30-50% relative difference) than backscatter 
with a minimum some 3-5% (30-60% relative) lower than the maximum. 

b. Konza 2 Site : Data from Ronza 2 provided an opportunity to examine 

the effect of sun zenith angle on bidirectional reflectance (Figs. 24 and 25). 
Reflectances in the visible wavebands 1 and 3 changed very little as the solar 
zenith angle decreased from 44 to 22°, particularly in the forward scattering 
direction. There was about a 1% change in reflectance (absolute terms) at 
nadir and at the backscatter angles. The effect of solar zenith angle on 
reflectance in band 4 was not clear. The highest reflectances for band 6 were 
measured at the lowest solar zenith angle. Canopy reflectances measured at 
the 34° solar zenith angle were consistently lower than reflectances measured 
at the 44° and 22° angles. Because reflectances can be influenced by changes 
in canopy and leaf structures and properties which could mask solar zenith 
effects measurements such as leaf water potential/ effective IAI and leaf 
angle distributions would provide useful information for helping to charac- 
terize changes in canopy reflectance over time. 

The effects of increasing IAI on bidirectional reflectance is demonstra- 
ted in Figs. 26-28. Three different leaf area coverages are illustrated with 
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LAI values of 0.5, 1.7 and 4.2. The contribution of soil reflectance to 
canopy reflectance decreased with increasing leaf area. This decrease is seen 
in the reflectance measurements. Visible and mid-IR reflectance decreased as 
leaf area increased while near-IR reflectance increased. Visible and mid-IR 
bidirectional reflectances did not differ between the areas with LAI of 1.7 
and those with LAI = 4.2 (Figs. 27 , 28). Areas with LAI of 1.7 apparently 
covered as much soil as those with LAI of 4.2 so that the reflection of 
radiation in these wavebands to the canopy reflectance was similar. 

An examination of all bidirectional reflectance data collected at Kbnza 
s ^ e ^ resulted in the following general conclusions: 

1) Reflectances in band 5 were generally higher than for band 4. This 
differs from the measurements made at Kbnza site 1 where reflectances 
were consistently higher for band 4. Changes in species composition 
or the influence of radiation from the soil may be a factor. 

2) Reflectance generally increased in the order of bands 1, 3, 2, 7, 6, 
4, 5. 

3) The minimum NIR reflectance occurred at 0 to 20° in the forward scat- 
tering direction with maximum reflectance at 50° in the backscatter 
direction. Forward scatter at 50° was generally 5-10% (10-20% rela- 
tive) lower than backscatter at the same angle. 

4) Minimum reflectance in the visible wavebands occurred in the 20 to 

30 forward scattering directions with the maximum in the 50° back- 

scatter direction. The minimum was generally 3-5% (35-49% relative) 
lower than the maximum. 

5) Maximum reflectance in channel 7 occurred in the 30 to 50° back- 
scatter direction and minimum in the 30 to 50° forward scatter 

direction. 



2. Albedo Calculations. 

Albedo is defined as the percentage of incident solar energy reflected 
away from a surface (lougeay, 1978; Rosenberg et al., 1983). It is an impor- 
tant quantity in the determination of the radiation balance of earth surfaces 
(Kriebel, 1979; Eaton and Dirmhim, 1979; Rosenberg et al., 1983). The impor- 
tance of albedo in the calculation of net radiation (Rn) is illustrated by 

equation ( 1 ) : 

Rn = (l-r)(Rsw+) + (Rlw+ - Rlwt) ^ 

where 

Rn = net radiation (W m"2) 

r = albedo (%) . . . _•>. 

Rswi = incoming shortwave radiation (W m ) 

Rlw4- = incoming longwave radiation (W m ) 

Rlwt = outgoing longwave radiation (W m" ) 

Albedoes may be different from one another and generally change with time 
(Kriebel, 1979). Albedo was expected to vary from place to place because the 
FIFE site presents variations in topography, management practices, vegetative 
conmunities, litter cover, soil exposure, and vegetation density. 

One of our ultimate goals is to calculate the net radiation of a surface 
employing data acquired by remote sensing instrumentation. The albedo must be 
known to quantify the shortwave balance (Eq. 1). Investigation into the uti- 
lity of employing remotely sensed data from the MR to calculate albedo was 
performed. These calculations of albedo utilized three components: 1) 

reflectance factors derived from the M© data, 2) the Walthall et al. (1985) 
bidirectional reflectance equation, and 3) the solar radiation curve of Moon 
(1940). The actual algorithm used to calculate the albedo is: 


7 

r = I Wi Rhi 

i=l 


(2) 


-A 5- 


where 


r = albedo 

Wi = a weighting factor 

Rhi = hemispherical reflectance for channel i as determined by the 
model of Walthall et al. (1985). 

A description of the reflectance factor calculations may be found in 
Gardner (1983) as elaborated upon by Bauer et al. (1981). (toon's (1940) solar 
radiation data provided solar irradiation in W m“ 2 jmT 1 for the region 0.295- 
2.14 ym for a surface at sea level with an optical air mass of two. The data 
of Moon were plotted and the NWR bandwidths were noted on the graph. 

To determine the albedo for a given surface using remotely sensed data, 
the bandwidths must account for all the radiation in the spectral curve 
constructed fran the Moon (1940) data. Therefore, each band of the (MR was 


extended 

to cover a little 

larger region of 

the spectrum than 

they actually 

measure 

(Table 4). 




Table 4. 

Weighting factors assigned to each MMR waveband based 
Moon (1940). 

on the data of 

Band 

Original 
Width ( ym) 

Extended 
Width (ym) 

Energy in 
Each Band 
(W m“ 2 ) 

Weighting 

Factor 

Wi 

1 

.45- .52 

.295- .520 

140 

0.1884 

2 

.52- .60 

.520- .615 

112 

0.1511 

3 

.63- .69 

.615- .700 

96 

0.1309 

4 

.76- .90 

.700-1.025 

228 

0.3072 

5 

1.15-1.30 

1.025-1.360 

106 

0.1435 

6 

1.55-1.75 

1.360-1.800 

49 

0.0665 

7 

2.08-2.35 

1.800-2.140* 

9 

0.0125 


*Moon's data showed no energy past 2.140 ym. 

Weighting factors are determined by dividing the energy in each band by 
the total energy under the curve (approximately 740 W m-2) (Moon, 1940). 


“ 46 “ 



Hemispherical reflectances (Phi) calculated for each channel for each sample 
plot are multiplied by their respective lighting factors (Wi) and then summed 
to provide the albedo for the sample plot. Albedos for a selected number of 
sample plots from the Konza 2 site were calculated and are shown in Table 5. 

Table 5. Albedo values calculated by equation 2 for several plots at the 
of f “Konza site* 


Plot # Time 


I 1026 24 

1 1126 21 

1 1247 24 

2 1034 25 

2 1119 22 

2 1255 25 

3 1044 20 

4 1053 21 

5 1101 22 

6 1108 21 

II 1340 25 


Estimated albedoes fall in the range of 20-25% which are reasonable 
albedo values for vegetative surfaces (Rosenberg et al., 1983). However, 
independent measurements of albedo for the sample plots were not made, there- 
fore, it is difficult to say that the determined values of albedo are correct. 

During the intensive field campaigns (IFCs) supporting measurements of 
net radiation, albedo, and canopy temperatures will be made in order to com- 
pare the albedo and net radiation values calculated from remotely sensed data. 
More recent and accurate solar radiation curves or models will be sought to 
improve albedo calculations. 

Objective (3) - Measurement of Sensible a nd Latent Heat Fluxes 

Measurements were made over a winter wheat crop (Triticum aestiv um L. , 
Colt) at Mead, Nebraska on selected days during the 1985 growing season. The 
incoming shortwave radiation (Fs) was measured with a pyrancmeter. Net radia- 
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tion (Rn) was measured using Swissteco net radiometers. Soil heat flux (S) 
was measured with heat flow sensors. Soil temperature was measured with ther- 
mocouples. Surface soil heat flux (S) was computed from soil heat flow sensor 
and soil thermocouple data employing a combination method. Mean wind speed 
profiles were measured with a set of three-cup anemometers. Profiles of mean 
air temperature and humidity were measured with ceramic wick psychrometers . 

Soil moisture was measured gravimetrically and with a Campbell Nuclear 
Pacific neutron probe. Crop height* leaf area index and dry matter were moni- 
tored. Water potential and stomatal resistance of flag leaves were also 
measured . 

Fluxes of sensible (H) and latent (LE) heat were computed by means of the 
Bowen-ratio technique. The data have been analyzed, tabulated and graphed. 
The following data sets were used to compare Cupid estimated fluxes with 
measured values. 


Variables 


a ) Meteorolog ical : 

Rs, Rn, S, H, LE, air temperature, 
vapor pressure deficit, wind speed, 
soil temperature 

b) Soil Moisture: 


c) Plant height, leaf width, height of 
most dense region, height of lowest 
leaves, leaf area index, dry matter: 


d) Leaf water potential and stcmatal 
resistance: 


Measurement Dates (1985) 

May 22, 23, 25, 28, 29, 30; 
June 7, 8, 13, 14 and 20 


May 21, 23, 28; June 3, 6 
8, 10, 12, 17, 18 and 24 

April 20, 24; May 6, 14, 

17, 20, 22, 23, 24, 28, 29, 
30; June 3, 5, 7, 8, 10, 12 
13, 14, 17, 18, 21 and 24 

May 20, 21, 22, 23, 24, 28, 
29 , 30 , 31; June 3, 5, 6, 7, 
8 and 10 


Fluxes of water vapor and sensible heat were calculated using the eddy 
correlation technique for a tall grass prairie location within the FIFE site 
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near Manhattan, Kansas on a few selected days during July-August, 1986. Ver- 
tical velocity fluctuations were measured with a sonic anemometer, air tem- 
perature and hunidity fluctuations were measured with fine wire thermocouples 
and a Iyman-alpha hygrometer. These measurements were made 2.25 m above the 
ground. Supporting data on net radiation, mean horizontal wind speed pro- 
files, soil heat flux, mean air temperature and humidity, and soil temperature 
were also obtained. Micrcmeteorological data were collected when the wind was 
south to east, resulting in approximately 200 m of upwind fetch over reason- 
ably flat terrain, volumetric soil water was determined gravimetrically in 

the upper 0.6 m soil layer. 

Signals from the eddy correlation sensors were recorded on an IBM PCAT 
microcomputer using an analog to digital converter. These signals were low- 
pass filtered with 8-pole Butterwerth active filters (12.5 Hz cutoff frequen- 
cy) and recorded at 25 Hz. The slaw-response signals were recorded on an IBM 
PCXT microcomputer . Data were averaged for the first 45 min of each solar 
hour. The last 15 minutes were reserved for sensor calibration and checks. 

Energy Fluxes 

The energy balance over the prairie vegetation can be approximated by 1 : 

Rn + IE + H + S- 0 

Values of Rn and S were treasured with slow-response sensors (net radiometers 
and heat flow transducers). On the other hand, IE and H were measured with 
fast-response eddy correlation sensors. Ideally, -(LE+H) should be balanced 
by (Rn+S). The sum of latent and sensible heat fluxes l- (LE+H)] measured by 


fluxes directed away from the surface are negative. 
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the eddy correlation technique generally agreed well with the sum of the net 
radiation and soil heat flux (Rn+S) , measured independently (Fig. 29). 

Diurnal patterns of the energy balance components over the tall grass 
prairie measured on 30 July 1986 are shown in Fig. 30a. The vegetation was at 
the heading stage and the leaf area index (LAI) was about 2.0. 

Except for a few early morning clouds, 30 July was clear with high solar 
radiation. Midday air temperature and vapor pressure deficit on this day 
ranged from 33 to 37°C and from 3.3 to 4.2 kPa, respectively (Fig. 31). 
Fluxes H, LE and S followed diurnal patterns similar to that of Rn. Net 
radiation reached its maximum of about 570 W m -2 at solar noon and decreased 
in the afternoon. Sensible and soil heat fluxes reached peak magnitudes of 
about 130 W m“ 2 and 90 W m -2 at solar nocan, respectively. Both fluxes were 
directed away from the surface throughout the day. The latent heat flux 
reached a peak magnitude of 410 W m" 2 around 1100 hrs and then decreased later 
in the day. 

The diurnal patterns of the energy balance components on a partly cloudy 
day (6 August 1986), at the same stage of growth, are shown in Fig. 30b. This 
day was mostly to partly cloudy in the morning hours and clear in the after- 
noon. Midday air temperature and vapor pressure deficit ranged from 27 to 
30°C and from 018 to 1.3 kPa, respectively. Values of Rn varied, depending on 
the degree of cloud cover, from 460 to 580 W nr 2 during the midday. Diurnal 
patterns of LE, H and S followed that of Rn quite closely. The latent heat 
flux was reduced due to lower Rn and atmospheric evaporative demand (lower air 
temperature, vapor pressure deficit and wind speed) on this day as compared to 
those on 30 July (see Fig. 31). A greater proportion of Rn was partitioned 
into sensible heat. Midday magnitude of H varied from 140 to 190 W m“ 2 . Mag- 
nitude of LE during the midday ranged between 210 to 340 W m“ 2 . Midday S was 
about 80 W m -2 * 
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Fig. 29. Balance of energy budget terms. LE and H measured by eddy 
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Profiles of volumetric soil water content in the upper 0.6 in, Treasured 
during 29 July - 6 August, are shown in Fig. 32. These data show that there 
was more water available in the soil on 6 August than on 30 July; therefore, 
the greater partitioning of Rn into H was not likely due to increased water 
stress at this later date. 


Objective (4) - Modeling Radiation, Sensible and Latent Heat Fluxes 

A. Estimated LAI and Mean Leaf Angles : 

Experimental and theoretical studies have demonstrated that the structure 
of a vegetative canopy can have a significant effect on its reflective pro- 
perties in the visible and NIR portions of the spectrum (Pinter 1985, cited by 
Jackson and Pinter 1986). Canopy architecture information can provide con- 
siderable aid in the interpretation and prediction of remotely sensed data. 
Direct measurements of canopy structure can be rather time consuming. The 
nature of the research to be conducted during FIFE will require a relatively 
rapid means of obtaining canopy structural measurements . 

An indirect non-destructive measure of LAI and leaf angle distribution 
was introduced by Norman et al. (1983). The method employs the relationship 
between the extinction of radiation beneath a canopy and the canopy structure. 
The relationship described mathematically is: 


T = e ~k ( 0 , a) LAI 


(4) 


T represents the transmission of radiation for a canopy, 0 is the sun or view 
zenith angle and k is the extinction coefficient given by: 


ir/2 

k(0, a) = / g(a) COS6 da (5) 

o 

g(a) is the fraction of leaves with inclination angle a to the horizontal, 6 
is the angle between the leaf normal and the given sun direction, and 0 is 
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the zenith angle of the sun. The penetration of radiation is dependent on the 
number and distribution of gaps in the canopy. Thus, canopy light penetration 
measurements contain a wealth of information on the vegetative structure. 

Using the relationship in Eq. (4) and given values of T and e, values of 
k and LAI can be inferred (Norman et al., 1979; Lang et al., 1985). Recently, 
Norman and Campbell (1987) proposed an inversion routine describing the canopy 
leaf angle distribution as an elliptical distribution. The ellipsoid can vary 
in horizontal and vertical axes describing most leaf angle distributions. The 
extinction coefficient, k, is defined as a function of LAI and x (the ratio of 
vertical to horizontal axes of the ellipse) and then solve for x and LAI such 
that: 


f - l (InTi + ki IAI) 2 is minimum. (6) 

The bisection method is used to find the solution where the RMS error is 
minimized. 

Canopy structure information has been inferred using photographic tech- 
niques and sunfleck measurements below a canopy with success in turf grass, 
com and soybean canopies (Lang et al., 1985; Perry, 1985; Kopeck, 1987). The 
view zenith angle is varied in the photographic method while the source zenith 
angle is varied in the sunfleck method. 

Cone sensors have been developed which sense diffuse light entering at a 
variety of angles (Norman et al., 1983). T\*o cone sensors were used in this 
research. One was a multi-cone sensor which measures diffuse light at angles 
12°, 33°, 49°, 65° and 81° (Norman et al., 1983). The second was a prototype 
sensor developed by LI-COR vhich measures diffuse light at 7°, 25°, 40° and 
58°. 

Tten measurements above and below the canopy were made in each plot. The 
LI-COR sensor required correction of raw data for corruption of light detected 
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by one or nore sensors. The ratio of corrected measurements below the canopy 
with neasurements above represent the fraction of diffuse light penetrating 
the canopy. The 10 ratios from each plot were averaged to give a mean gap 
fraction at each angle. These values were used as input values for the inver- 
sion routine. 

Part of the 1986 field experiment included testing the suitability of the 
multi-cone sensor in a grassland environment. LAI estimates inferred using 
the multi-cone sensor data were compared to IAI estimates inferred from the 
LI-GQR oone data. Multi-cone IAI estimates were generally Lower than LI-COR 
cone IAI estimates (Fig. 33). The lower values of the multi-cone IAIs are 
due, in part, to the height of the oone (approximately 0.1 m) so that the leaf 
area below the height of the cones was not viewed. Canopy heights were as low 
as 0.30 m in some plots. A shorter oone like the LI-COR instrument is 
desirable for grassland applications. Also the multi-cone sensor is more sen- 
sitive in the green portion of the spectrum, thereby sensing more scattered 
radiation. The canopy appears to be brighter at greater zenith view angles 
where it should be darker; this results in higher calculated values of gap 
fraction at these lower angles. The distribution of gap fractions with cone 
angle were exemplary of canopies with lower leaf area and more horizontally 
inclined leaves than actually existed. In most cases, the multi-cone sensor 
data resulted in lower IAI and lower mean leaf angle estimates than those 
derived from the LI— OOR sensor. Estimates of IAI, as well as additional 
canopy structural measurements, are given in Table 6. 

Inferred IAI was plotted against dry leaf weight and canopy height (Figs. 
34, 35). The relationship is described with an exponential equation from 
Landsberg (1977) of the forms 

IAI(s) = A(1 - exp B * X+C *) 
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Table 6. Canopy structure characteristics, FIFE pilot study. July, August, 
1986. 


Plot # 

Location 

Dry 

Weight 

(g) 

Height 

(an) 

Inferred 

LAI 

Mean 

Leaf 

Angle 

1 

KONZA 2 

64.0 

35 

2.4 

62.5 

2 

KONZA 2 

81.4 

35 

2.6 

62.5 

3 

KONZA 2 

60.1 

30 

1.1 

74.9 

4 

KONZA 2 

64.5 

30 

1.3 

69.3 

5 

KONZA 2 

92.1 

40 

2.4 

64.1 

6 

KONZA 2 

71.1 

35 

3.1 

59.4 

7 

KONZA 2 

149.8 

65 

5.1 

50.2 

8 

KONZA 2 

131.8 

45 

3.4 

53.3 

9 

KONZA 2 

113.6 

50 

4.4 

57.9 

10 

KONZA 2 

93.7 

40-45 

3.7 

55.8 

11 

KONZA 2 

104.8 

30-35 

3.4 

64.1 

21 

KONZA 1 

106.7 

35 

2.6 

46.0 

22 

KONZA 1 

77.7 

35 

2.0 

50.9 

23 

KONZA 1 

160.5 

70 

3.7 

56.7 

24 

KONZA 1 

195.9 

70 

4.5 

55.0 

25 

KONZA 1 

200.4 

65 

2.9 

62.5 

26 

KONZA 1 

114.4 

65 

3.5 

62.5 

27 

KONZA 1 

129.8 

60 

3.4 

54.1 

28 

KONZA 1 

119.0 

70-75 

3.6* 


29 

KONZA 1 

93.7 

50 

2.9* 


30 

KONZA 1 

137.1 

60 

4.8* 


31 

KONZA 1 

112.5 

45 

3.5* 


32 

KONZA 1 

57.3 

30 

1.4* 



* Approximate LAI; values extrapolated frcm curve 
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Figure 35. 



where x is the leaf or total biomass dryweight in grams or canopy height in 
cm. Such relationships can be useful in estimating LAI from the dry biomass 
or, conversely, estimating biomass from estimates of LAI which can be made 
quite rapidly using sensors such as the one developed by LI -COR. 

B. Model Simulations of Wheat Canopy Fluxes : 

Coordinated measurements of biological characteristics, hourly climatic 
data and fluxes made over a wheat canopy at Mead, Nebraska in May and June, 
1985 and over a site on the FIFE pilot study area periodically from July 29 
through August 6, 1986 provided opportunities to test the performance of the 
model, Cupid (Norman, 1982) for agronomic and natural grassland communities. 

Cupid is a comprehensive plant-environment model, including many environ- 
mental and plant factors. Ihe model requires a number of input values which 
make the model versatile. Input information includes hourly measurements of 
solar radiation, air temperature, relative humidity, wind speed, precipita- 
tion and soil temperature. Hourly inputs were available from the automated 
weather stations located at the study areas. Canopy input information on 
soils and vegetation makes the model independent of surface cover. Soil 
inputs include bulk density, soil texture, volumetric soil moisture at seme 
reference depths, soil emissivity and soil reflection properties in the 
visible, near— infrared (NIR) and thermal portions of the spectrum. Vegetation 
input includes leaf area index, canopy height, mean leaf angle, leaf reflec- 
tance and transmittance (in the visible, near-NIR and thermal), and vegetation 
emissivity. Stcmatal characteristics for the wheat canopy were taken from 

Denmead and Millar (1976) while those for the FIFE site were determined from 
field measurements. 
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C. Vfoeat Canopy Fluxes ; 


We examined measured and predicted fluxes of Rn, I£» H and S over a wheat 
canopy for four days in 1985: May 22 (Day 142), May 30 (Day 150), June 7 (Day 
158) and June 13 (Day 164). These represent four days of contrasting weather 

and canopy conditions. 

Soils at Mead are a Sharpsburg silty clay loam. Wheat row spacing was 
0.18 m with a plant spacing of 0.03 m. Leaf area index of the canopy 
decreased with time. The LAI for the four days were 4.5, 3.0, 1.2 and 0.3, 
respectively. Radiation, wind, temperature and vapor pressure input were 
available from the Mead automated weather station, providing hourly input over 
a 24-hour period (Figs. 36-39). No precipitation was recorded on any of these 

days. 

There was very good agreement between measured and simulated Rn, l£ and H 
fluxes for May 22 and May 30 (Figs. 40-43), clear and partly cloudy days with 
high but decreasing LAI. Measured and simulated fluxes of Rn for June 7 and 
June 13 were also in good agreement although simulated IE and H fluxes were in 
poor agreement (Figs. 44-47). The poor agreement probably results from the 
senescing stages of wheat and a significant decrease in green leaf area. 
Although LAI had decreased considerably by day 164, wheat heads probably 
contributed more to the biomass than was represented by the model. Charac- 
terizing the wheat heads in the model should improve simulations. 

D. Natural Grassland Canopy Fluxes : 

Measured and predicted fluxes of Rn, IE, H and S for two days: July 30 
(Day 211) and August 3 (Day 215) were examined. These represent days of 
contrasting weather conditions. LAI was about 2.9 on both dates. July 30 was 
mostly sunny with temperatures as high as 40°C (Fig. 48). August 3 was cloudy 
near solar noon with lower temperatures than July 30 (maximum - 29.2°C) (Fig. 
49). Winds of 6 m/s were recorded during the early morning hours. 
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WHEAT — MEAD, NE 5x22x85 (DAY M 2 ) 



HOUR CCST) 


Figure 40. Caparison of measured and simulated Rn and IZ fluxes on May 22. 1985 (Day 142) 
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Figure 41. Comparison of measured and simulated H fluxes on May 22, 1985 (Day 142) 
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Figure 42. Comparison of 


measured and simulated Rn and 


LE fluxes on May 30, 1985 


(Day 150). 



Figure 44. Comparison of measured and simulated Rn and LE fluxes on June 7, 1985 (Day 158 
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Figure 45. Comparison of measured and simulated H fluxes on June 7, 1985 (Day 158). 



FLUX (W/m 


WHEAT — MEAD, NE 6/ 13/85 (DAY 164) 



0 G 12 18 24 

HOUR ( CST ) 

Figure 46. Comparison of measured and simulated Rn and IE fluxes on June 13, 1985 (Day 164). 
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Figure 47. Comparison of measured and simulated H fluxes on June 13, 1985 (Day 164). 
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Figure 48. Hourly values of solar radiation, wind speed and 
air temperature on July 30, 1986 (Day 211). 
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Figure 49. Hourly values of solar radiation, wind speed and 
air temperature on August 3, 1986 (Day 215). 
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Simulations of Rn were in good agreement with measured values for all 
days (Figs. 50 and 51). LE was overestimated for each case , although trends 
between measured and simulated were similar although differing in magnitude. 
H was consistently underestimated and was in poor agreement on August 3 (Figs. 
52 and 53). Simulated values of S were in good agreement (Figs. 54 and 55). 

CONCLUSIONS 

Estimates of Rn, LE and H made in 1985 with the model Cupid agreed well 
with the estimates of these fluxes made with micrometeorological techniques 
over wheat, especially in May before plant senescence had begun. The agree- 
ment was less good in June as LAI values decreased significantly. The perfor- 
mance of the Cupid model for estimating LE, Rn and H fluxes in 1986 fran 
prairie vegetation was poorer than that for the more homogeneous wheat canopy 
in 1985 even at comparable LAI values. Modifications to Cupid are needed to 
account for the nonhomogeneous, randomly-spaced clumps of prairie vegetation 
as well as a need to incorporate the influence of leaf litter on the fluxes of 
LE, H and S. 
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Figure 50. Measured and simulated values of net radiation 
(Rn) and latent heat flux (IE) for July 30, 1986 
(Day 211). 
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Figure 51 . Measured and simulated values of net radiation 
(Rn) and latent heat flux (LE) for August 3, 
1986 (Day 215). 


- 80 - 






Measured S 


'B i mu 1 a. ted S 



0 4 8 12 16 20 24 

HOUR OF DHY 


Figure 54. Measured and simulated values of soil heat flux 
(S) for July 30, 1986 (Day 211). 
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Figure 55. Measured and simulated values of soil heat flux 
(S) for August 3, 1986 (Day 215). 
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